Controlling Molecular Scattering by Laser-Induced Field-Free Alignment 
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We consider deflection of polarizable molecules by inhomogeneous optical fields, and analyze the 
role of molecular orientation and rotation in the scattering process. It is shown that molecular rota- 
tion induces spectacular rainbow-like features in the distribution of the scattering angle. Moreover, 
by preshaping molecular angular distribution with the help of short and strong femtosecond laser 
pulses, one may efficiently control the scattering process, manipulate the average deflection angle 
and its distribution, and reduce substantially the angular dispersion of the deflected molecules. We 
provide quantum and classical treatment of the deflection process. The effects of strong deflecting 
field on the scattering of rotating molecules are considered by the means of the adiabatic invariants 
formalism. This new control scheme opens new ways for many applications involving molecular 
focusing, guiding and trapping by optical and static fields. 

PACS numbers: 33.80.-b, 37.10.Vz, 42.65.Re, 37.20.+j 



I. INTRODUCTION 

Optical dipole forces acting on molecules in nonreso- 
nant laser fields is a hot subject of many recent exper- 
imental studies By controlling molecular transla- 
tional degrees of freedom with laser fields [a-llO]. novel 
elements of molecular optics can be realized, includ- 
ing molecular lens [3, [|| and molecular prism Q . The 
mechanism of molecular interaction with a nonuniform 
laser field is rather clear: the field induces molecular po- 
larization, interacts with it, and deflects the molecules 
along the intensity gradient. As most molecules have 
anisotropic polarizability, the deflecting force depends on 
the molecular orientation with respect to the deflecting 
field. Previous studies on optical molecular deflection 
have mostly considered randomly oriented molecules, for 
which the deflection angle is somehow dispersed around 
the mean value determined by the orientation-averaged 
polarizability. The latter becomes intensity-dependent 
for strong enough fields due to the field-induced modifi- 
cation of the molecular angular motion [5|, This adds 
a new ingredient for controlling molecular trajectories 0- 
■8], which is important, but somehow limited because of 
using the same fields for the deflection process and ori- 
entation control. 

In this work, we show that the deflection process can be 
significantly affected and controlled by preshaping molec- 
ular angular distribution before the molecules enter the 
interaction zone. This can be done with the help of nu- 
merous recent techniques for laser molecular alignment, 
which use single or multiple short laser pulses (transform- 
limited, or shaped) to align molecular axes along certain 
directions. Short laser pulses excite rotational wavepack- 
ets, which results in a considerable transient molecular 
alignment after the laser pulse is over, i.e. at field-free 
conditions (for reviews on field-free alignment, see, e.g. 
HI!])- Field- free alignment was observed both for 
small diatomic molecules as well as for more complex 
molecules, for which full three-dimensional control was 
realized 0-Qjl. 

We demonstrate that the average scattering angle of 



the deflected molecules and its distribution may be dra- 
matically modified by a proper field-free prealignment. 
By separating the processes of the angular shaping and 
the actual deflection, one gets a flexible tool for tailoring 
molecular motion in inhomogeneous optical and static 
fields. 

The main principles of this new approach were briefly 
introduced in our recent Letter [17] . Here we present 
a much more elaborated analysis of the control mecha- 
nisms, including also a detailed comparison between the 
quantum and classical aspects of the problem, and dis- 
cussion of the strong field effects in molecular scattering. 

In Sec. [IT] we present the deflection scheme, as well 
as heuristic classical discussion on the anticipated role 
of molecular rotation on the deflection process (both for 
thermal and prealigned molecules) . In Sec IIIII we verify 
these predictions by the means of the quantum treatment 
of the problem in the limit of the relatively weak deflect- 
ing field (that does not disturb significantly the rotational 
motion). The strength of the prealigning field is not re- 
stricted here. Full classical treatment of the molecular 
deflection at such conditions (including thermal effects) 
is given in Sec. IIV[ where we find a good correspondence 
between the classical and quantum calculations. Moti- 
vated by this agreement, we provide in Sec. [V] a full 
classical analysis of the molecular scattering by strong 
deflecting field using the adiabatic invariants formalism. 
Finally we summarize our results in Sec. [VI] 



II. DEFLECTION OF FIELD-FREE ALIGNED 
MOLECULES 

Although our arguments are rather general, we follow 
for certainty a deflection scheme that brings to mind the 
experiment by Stapelfeldt et al [l[ who used a strong 
IR laser to deflect a CS2 molecular beam, and then ad- 
dressed a portion of the deflected molecules (at a prese- 
lected place and time) by an additional short and narrow 
ionizing pulse. Consider deflection (in z direction) of a 
linear molecule moving in x direction with velocity v x 
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and interacting with a focused nonresonant laser beam 
that propagates along the y axis (Fig. [!}. 




FIG. 1: The deflection scheme. A polarized (in the z direc- 
tion) laser field propagates toward the plane of the paper (y 
direction). The linear molecules, initially moving along the 
x direction (with velocity v x ), axe deflected by the potential 
gradient (deflection velocity v z ). 

The spatial profile of the laser electric field in the xz- 
plane is: 

E = E exp[~(x 2 + z 2 )/uj 2 } exp[-21n2i 2 /T 2 ]. (1) 

The interaction potential of a linear molecule in the laser 
field is given by: 

U = -^E 2 (a,, cos 2 9 + a ± sin 2 9) , (2) 

where E is defined in Eq. [TJ and ay and a± are the com- 
ponents of the molecular polarizability along the molec- 
ular axis, and perpendicular to it, respectively. Here 9 is 
the angle between the electric field polarization direction 
(along the laboratory z axis) and the molecular axis. A 
molecule initially moving along the x direction will ac- 
quire a velocity component v z along z-direction. We con- 
sider the perturbation regime corresponding to a small 
deflection angle, 7 w V z /v x . We substitute x = v x t, and 
consider z as a fixed impact parameter. The deflection 
velocity is given by: 




Here M is the mass of the molecules, and F z is the de- 
flecting force. The time-dependence of the force F z (and 
potential U) in Eq.Q comes from three sources: pulse 
envelope, projectile motion of the molecule through the 
laser focal area, and time variation of the angle 9 due to 
molecular rotation. For simplicity, we start with the case 
of the relatively weak deflecting field that does not affect 
significantly the rotational motion. Such approximation 
is justified, say for CS2 molecules with the rotational 
temperature T — 5K, which are subject to the deflect- 
ing field of 3 • 10 9 W/cm 2 . The corresponding alignment 
potential U w — \ (an — a±) Eq « 0.04 meV is an or- 
der of magnitude smaller than the thermal energy ksT, 
where fc^ is Boltzmann's constant. This assumption is 
even more valid if the molecules were additionally sub- 
ject to the aligning pulses prior to deflection. The case of 
a strong deflecting field will be considered later in SeclVl 

Since the rotational time scale is the shortest one in 
the problem, we average the force over the fast rotation, 
and arrive at the following expression for the deflection 
angle, 7 = v z /v x : 

7 = 70 [a\\A + ot±(l- A)] /a (4) 

Here a = l/3Q!i| + 2/3a± is the orientation- averaged 

molecular polarizability, and A = cos 2 9 denotes the 
time-averaged value of cos 2 9. This quantity depends on 
the relative orientation of the vector of angular momen- 
tum and the polarization of the deflecting field. It is 
different for different molecules of the incident ensem- 
ble, which leads to the randomization of the deflection 
process. The constant 70 presents the average deflection 
angle for an isotropic molecular ensemble: 

aE 2 f-4:z\ 
70 ~ AMvl \~ ) 

R ( 2w 2 ln2V 1/2 / 2z 2 \ , . 

We provide below some heuristic classical arguments on 
the anticipated statistical properties of A and 7 (both for 
thermal and prealigned molecules). 

Consider a linear molecule that rotates freely in a 
plane that is perpendicular to the vector J of the 
angular momentum (see Fig.©). 

The projection of the molecular axis on the vertical 
z-direction is given by: 

cos 9(t) = cos(wi) sin^j, (6) 

where 9j is the angle between J and z-axis, and u is the 
angular frequency of molecular rotation. Averaging over 
time, one arrives at: 

A = ~c~^ r 9 = \ sin 2 9j. (7) 
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FIG. 2: A molecule rotates with the angular momentum J 
forming an angle 8j with the laboratory z axis. 



In a thermal ensemble, vector J is randomly oriented 
in space, with isotropic angular distribution density 
l/2sin(0j). The mean value of the deflection angle is 
(7) = Yo. Eq.([7]) allows us to obtain the distribution 
function, f(A) for A (and the related deflection angle) 
from the known isotropic distribution for 9j. Since the 
inverse function Oj(A) is multivalued, one obtains 
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where we summed over the two branches of 6j(A). 
This formula predicts an unimodal rainbow singularity 
in the distribution of the scattering angles at the 
maximal value 7 = 70 (en 11 + a±)/2a (for A — 1/2), and 
a flat step near the minimal one 7 = joa±/a (for .4 = 0). 

Assume now that the molecules are prealigned before 
entering the deflection zone by a strong and short laser 
pulse that is polarized perpendicular to the polarization 
direction of the deflecting field (e.g., in x-direction). Such 
a pulse forces the molecules to rotate preferentially in the 
planes containing the x-axis. As a result, the vector J of 
the angular momentum is confined to the yz-plane, and 
angle 0j becomes uniformly distributed in the interval 
[0,7r] with probability density I/71". The corresponding 
probability distribution for A takes the form 



f(A) 



V2 



1 



7T y/A(l ~ 2A) 



(9) 



In contrast to Eq.([5]). Eq.© suggests a bimodal rainbow 
in the distribution of deflection angles, with singularities 
both at the minimal and the maximal angles. Finally, we 
proceed to the most interesting case when the molecules 
are prealigned by a short strong laser pulse that is polar- 
ized parallel to the direction of the deflecting field. After 
excitation by such a pulse, the vector of the angular mo- 
mentum of the molecules is preferentially confined to the 
xy-plane, and the angle 8j takes a well defined value of 
9 j k, 7r/2 which corresponds to A = 1/2. In this way, 
the molecules experience the maximally possible time- 
averaged deflecting force which is the same for all the 
particles of the ensemble. As a result, the dispersion of 
the scattering angles is reduced dramatically. The distri- 
bution of the deflection angle 7 transforms to a narrow 



peak (asymptotically - a <5-function) near the maximal 
value, 7 = 70 (cum + a±)/2a. 



III. QUANTUM TREATMENT 

For a more quantitative treatment, involving analysis 
of the relative role of the quantum and thermal effects on 
one hand, and the strength of the prealigning pulses on 
the other hand, we consider quantum-mechanically the 
deflection of a linear molecule described by the Hamilto- 



n = J 2 7(27). 



(10) 



Here J is operator of angular momentum, and 7 is the 
moment of inertia, which is related to the molecular ro- 
tational constant, B — h/(4irlc) (c is speed of light). 
Assuming again that the deflecting field is too weak to 
modify molecular alignment, we consider scattering in 
different | J, m) states independently. The deflection an- 
gle is given by Eq.Q, in which A is replaced by 



4j, m = (J,m\ cos 2 



1 2 7(7 + 1) - 3m 
1 ' TO)_ 3 + 3(2J + 3)(2J 



1) 



(11) 

In the quantum case, the continuous distribution of the 
angles 7 is replaced by a set of discrete lines, each of 
them weighted by the population of the state | J, m) . 
Fig. [3] shows the distribution of Aj, m in the thermal 
case for various values of the dimensionless parameter 
Jt = \J ksT/ (hBc) that represents the typical " thermal" 
value of J (for Jt > 1). For CS2 molecules, the values 
of J T = 5, 15 correspond to T = 3.9K and T = 35K, 
respectively. 

The distribution of discrete values of Aj im demon- 
strates a non-trivial pattern. In particular, the values 
exceeding the classical limit 0.5 correspond to the states 
1 7, m = 0) (see Eq. fTTj) ). and they rapidly approach that 
limit as J grows. After the coarse-grained averaging, 
however, the distribution shows the expected unimodal 
rainbow feature (see Eq.©) for large enough Jt- 

If the molecules are subject to a strong femtosecond 
prealigning pulse, the corresponding interaction poten- 
tial is given by Eq. ([2]), in which E is replaced by the 
envelope e of the femtosecond pulse. If the pulse is short 
compared to the typical periods of molecular rotation, it 
may be considered as a delta-pulse. In the impulsive ap- 
proximation, one obtains the following relation between 
the angular wavefunction before and after the pulse ap- 
plied at t = (see e.g. [18j . and references therein): 

#(t = 0+) = exp(iFcos 2 0)*(t = 0"), (12) 

where the kick strength, P is given by: 



P= (1/471) • (ck|| -ol) / e z {t)dt. (13) 
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FIG. 3: Quantum distribution of Aj, m in the thermal case. 
Panels (a) and (b) correspond to Jt = 5 and Jt = 15, respec- 
tively. Histogram in panel (c) shows a coarse-grained version 
of the distribution in panel (b). 



Here we assumed the vertical polarization (along z-axis) 
of the pulse. Physically, the dimensionless kick strength, 
P equals to the typical amount of angular momentum 
(in the units of ft) supplied by the pulse to the molecule. 
For example, in the case of CS2 molecules, the values of 
P = 5,25 correspond to the excitation by 0.5ps (FWHM) 
laser pulses with the maximal intensity of 4.6- 10 n Wy cm 2 
and 2.3 • 10 12 W/cm 2 , respectively. For the vertical po- 
larization of the laser field, m is a conserved quantum 
number. This allows us to consider the excitation of the 
states with different initial m values separately. In or- 
der to find ^(t = + ) for any initial state, we introduce 
an artificial parameter £ that will be assigned the value 
£ = 1 at the end of the calculations, and define 

4\ =cxp[(iPcos 2 0)£] *(t = - ) = ^c/(£)|J,m). 

(14) 

By differentiating both sides of Eq. (TT4"|) with respect to 
£, we obtain the following set of differential equations for 
the coefficients cj\ 



cji = iP^2cj(J', m\ cos 2 0\J, to), 



(15) 



where c = dc/d^. The diagonal matrix elements in 
Eq. (|15p are given by Eq. (fTTj) . the off-diagonal ones can 
be found using recurrence relations for the spherical har- 
monics [13]. Since * 5=0 = #(t = 0~) and * 4= i = ^(t = 
+ ) (sec Eq. (fT4|) ). we solve numerically this set of equa- 
tions from £ = to £ = 1, and find ^ (t — Q + ) . In order to 
consider the effect of the field-free alignment at thermal 
conditions, we repeated this procedure for every initial 



I Jo, mo) state. To find the modified population of the 
I J, to) states, the corresponding contributions from dif- 
ferent initial states were summed together weighted with 
the Boltzmann's statistical factors: 



f(Aj, mo ) = Yl 

Jo, J 



cxp{-E Jo /K B T) 



x \cj\ 2 5(Aj,m ,Aj 7mo ), 



(16) 



where cj are the coefficients (from Eq. [T5)) of the wave 
packet that was excited from the initial state | Jo, too); 5 is 
the Kronecker delta, and Q ro t is the rotational partition 
function. It worth mentioning that different combina- 
tions of J and to may correspond to the same value of 
Aj >m , which necessitates the presence of the Kronecker 
delta in Eq. (fFo)) . For symmetric molecules, statistical 
spin factor should be taken into account. For example, 
for CS2 molecules in the ground electronic and vibra- 
tional state, only even J values are allowed due to the 
permutation symmetry for the exchange of two Bosonic 
Sulfur atoms (that have nuclear spin 0). 

In the case of an aligning pulse in the x direction, the 
operator in Eq. fT!?)) becomes: 



*(i = 0+) =exp(iPcos 2 0sin 2 6»)*(t = 0" 



(17) 



and a similar procedure as described above is used to 
find the deflection distribution. One should pay attention 
that to is no longer a conserved quantum number for a 
pulse kicking in the x direction. 




FIG. 4: Distribution of A.j, m for molecules prealigned with 
the help of a short laser pulse polarized in the x direction. 
The left column (a-b) presents directly the Aj, m values, while 
the right column (c-d) shows the corresponding coarse-grained 
histograms (as in Fig. [3J:). Panels (a) and (c) are calculated 
for Jt — 5 and P = 5; (b) and (d) are for Jt = 5 and P — 25. 

Using this technique, we considered deflection of ini- 
tially thermal molecules that were prealigned with the 
help of short pulses polarized in x and z directions (Figs. 
2] and respectively) . In the case of the alignment per- 
pendicular to the deflecting field (Fig. 0|, the coarse- 
grained distribution of Aj, m (and that of the deflection 
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which leads to 



FIG. 5: Distribution of Aj, m for molecules prealigned in the 
z direction. The left column (a-b) presents directly the Aj,m 
values, while the right column (c-d) shows the corresponding 
coarse-grained histograms. Panels (a) and (c) are calculated 
for Jt = 5 and P — 5; (b) and (d) are for Jt = 5 and P = 25. 



angle) exhibits the bimodal rainbow shape, Eq.® for 
strong enough kicks (P 3> 1 and P ^> Jt). Finally, 
and most importantly, prealignment in the direction par- 
allel to the deflecting field allows for almost complete re- 
moval of the rotational broadening. A considerable nar- 
rowing of the distribution can be seen when comparing 
Fig. [3^ and Figs. [5p and[5jl. The conditions required 
for the considerable narrowing shown at Fig. [5ji corre- 
spond to the maximal degree of field-free pre-alighment 
(cos 2 9)max = 0.7. This can be readily achieved with the 
current experimental technology, even at room tempera- 
ture El 



IV. CLASSICAL TREATMENT: WEAK 
DEFLECTING FIELD 

Consider a classical rigid rotor (linear molecule) de- 
scribed by a Lagrangian 



b 2 sin 2 ( 



(18) 



where 9 and cj> are Euler angles, and I is the moment of 
inertia. The canonical momentum for the 4> angle 



Pa, = sin 2 



(19) 



is a constant of motion as is a cyclic coordinate. The 
canonical momentum Pq is given by 

Pe = 16 (20) 

The Euler-Lagrange equation for the 9 variable is 



d dL dL 
dt'dd ~d9 



0, 



(21) 



-P 2 cos6> 
~Pl^9 



(22) 



When considering a thermal ensemble of molecules, it is 
convenient to switch to dimensionless variables, in which 
the canonical momenta are measured in the units of pth = 
Iuhhi with Lu t h = y/kBT/I, where T is the temperature 
PH. By setting P£ = P^/p th , P' e = Pel Pth, and f = 
cutht, one gets the following solution of Eq. 



cos(9(Q)+ut'), (23) 



where 



if (0) 



^ 2 (Q) ' 

sin 2 9(0) 



1/2 



(24) 



As in Sec IIII1 if at t = the molecules are subject to a 
femtosecond aligning pulse polarized in z-direction, the 
corresponding interaction potential is given by Eq. ([2]l. 
in which E is replaced by the envelope e of the femtosec- 
ond pulse. We assume again that such a pulse is short 
compared to the rotational period of the molecules, and 
consider it as a delta-pulse. The rotational dynamics of 
the laser-kicked molecules is then described by Eq. 
in which P' g (0) is replaced by 



fftO) -> P' e {Q) - P s 'sin(20(O)). 



(25) 



Here P' s = PK/ y/ksTI is properly normalized kick 
strength [l|] with P given by Eq. (fl«|. 

In the case of an aligning pulse in the x direction, both 
Pg{0) and PL(Q) are replaced by: 

P^(0) -> P e '(0) + P s ' cos 2 0(0) sin(26»(0)) 

P;(0) ^ P^(O)-F s 'sin 2 (0(O))sin(20(O)) (26) 

Averaging cos 2 9(f) over time, we obtain: 



A 



1 - 



P^(0) 



cos(26»(0)), (27) 



and the probability distribution of the time-averaged 
alignment factor A can be obtained by: 



f(A) = J J J J d9(0)d(j>(0)dP g '(0)dP^(0)S(A-cos 2 e) 
x F[0(O),0(O),P^(O),F;(O)] (28) 



where 



^=8^ 6XP 




(29) 
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is the thermal distribution function. 

The probability distribution of A in the thermal case 
is plotted in Fig. [5] Its shape is well described by Eq. [FJ 
and it is in good agreement with the quantum result of 
Fig. [3b. 




FIG. 6: Classical distribution of A for Jt = 15. No prealign- 
ment is assumed (P = 0). One can observe a rainbow- like 
feature at the right edge of the distribution, and a flat step 
at the left edge. 

Figs. [7] and [8J show the distribution of A value for 
molecules that were prealigned in the direction perpen- 
dicular and parallel to the deflection held, respectively. 
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FIG. 7: Classical distribution of A for (a) Jt = 0.5 and (b) 
Jt — 5 after the molecules were prealigned by a laser pulse 
(P = 25) in the x direction (i.e. perpendicular to the deflect- 
ing field). Figure (a) corresponds to the case P ^ Jt, and it 
is in good agreement with the analytical result, Eq. (9]). 

In the case of perpendicular prealignment by a suffi- 
ciently strong kick (P 3> Jt), the distribution shown in 
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FIG. 8: Narrow classical distribution of A results from pre- 
alignment by means of a laser pulse polarized parallel to the 
deflecting field (P = 25 and J T = 5). 



Fig. [7^ demonstrates the bimodal rainbow shape pre- 
dicted by Eq. © . Figure [7]d is similar to the corre- 
sponding quantum histogram in Fig. 0Ji. 

In the case of parallel prealignment, the predicted nar- 
row distribution is seen in Fig. [5J In what follows, we 
provide an asymptotic estimate of the width, AA of this 
distribution, and the mean value < A > of A in the limit 
of P/J T > 1. 

For strong enough kicks, Eq. (|27p shows that A ap- 
proaches the value of \, unless 0(0) is close to \. There- 
fore, we define 



5A = A-~. 



(30) 



that is different from zero only for small values of /? = 
0(0) - f . For 0(0) w f , Eqs. (25]), @7) and (00) yield: 



6 A: 



(p^(o) + 2/?p;r + ^(o)2 

The thermal averaging provides: 



(31) 



(SA) 



dP' 9 {Q)dP'^) 



/V^(0) + 2/3P s ') 2 + P'(0)2 



exp • 



{-\[pm 2 +pm 2 ]} 



(32) 



where Q ro t is the rotational partition function, and 
[— £, +£] is the interval of the f3 values for which the ap- 
proximation (I31[) is valid . We continue manipulating the 
expression by introducing 7 = 2/3 P' s : 



7 



and 



(SA) 



1 

4P 7 , 

-ap« (^(o)+7r+n(°) 2 

exp{-l[^(0) 2 + P;(0) 2 ]} 



(33) 



In the limit of P' s — > oo, the leading term in the asymp- 
totic expansion of (<L4) can be obtained by expanding 
the limits of the internal integration ±2P S '£ to ±oo (as 
the integrand vanishes for large values of P' s £,): 



(SA) 



4P' 



dpi{0)dP^0)\P'^)\ 



7 2 + l 



exp • 



{-i[^(0) 2 + P;(0) 2 ]} 



1 

Recalling Eq. (|30[) . we conclude that 



(34) 



(35) 



In order to estimate the width of the A distribution, 
we need to consider the dispersion, and accordingly the 
average value of A 2 . Following the same procedure as 
above, we define: 



S(A 2 ) = A 2 -\. 



(36) 



and find: 



5(A 2 ) 



P^(0) 4 + 2P^(0) 2 (P,'(0)+2/3P s ') 



(P^(0) + 2/?P0 2 + P'(0) 2 



By thermally averaging this, and taking only the leading 
term in the asymptotic expansion, we arrive at 



(A 2 ) 



16P' 



(39) 



The variance can be calculated from Eqs. (|35|) and ([39]) 
by using (A.4) 2 = (A 2 ) — (A) 2 . Recalling the relations 
P' s = Ph/y/k^TI and J T = y/k B T/(hBc), we have: 



(A) 
(AA) 2 



1 _ ^/tt J T 

2 ~ ~~8~ P~ 

Jt 

"32" ~p ' 



(40) 



The above asymptotic expressions for AA and (A) are 
plotted in Figures [9^l and [9]b, respectively (solid lines). 
The x points refer to the direct numerical calculations 
based on the distribution function given by Eg. (|2"51) . Al- 
though the asymptotic results (|4"0"j) are formally valid 
for Jt > 1, and P 3> Jt, they provide a good agree- 
ment with the exact numerical simulations already for 
P/Jt — 2. Moreover, our classical asymptotic estimate 
for the width of the distribution, AA coincides within the 
10% accuracy with the exact quantum result for P = 25 
and Jt = 5 presented above. 
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FIG. 9: Analytical asymptotic estimations of (a) A.4 and (b) 
(A) (solid lines) . The x points result from a direct numerical 
calculation using Eq. (|28[) . 



(6(A 2 )) 



1 

8P7 



6^(0)^(0)^(0)1 



l + 2 7 2 

7 ',,2 1 1W 



exp 



(7 2 + 1) 

{-^[^(Q) 2 + n(°) 2 ]} 



16P' 



Qrot 

n 
2' 



(38) 



V. CLASSICAL TREATMENT: STRONG 
DEFLECTING FIELD 

In the case of a strong deflecting field, the rotational 
motion of the molecules may be disturbed by the field. 
The Hamiltonian of a molecule in the vertically polarized 
optical electric field of a constant amplitude E is: 



H 



b 2 sin 2 I 



1 



E 2 (Aa cos 2 9 + a±) (41) 
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The conjugate momenta and Pg are given by Eqs. 
([H?| and ([20]) , where P^ is a constant of motion for the 
chosen polarization. It is convenient 21] to introduce a 
new variable 



u = cos 2 I 



(42) 



that satisfies the equation 



(1-«)i9- 



P,<, 



+ (1 — u)au 



The coefficients a,/3 in the polynomial g(u) are given by 
/^I^F+^V) (44) 

and 

ArvR 2 

(45) 



2/ 

Equation (|4"3")l can be immediately solved by separation 
of variables 



dt = 



du 



(46) 



In the case of a free rotation (a — 0), there are only 
two roots to the polynomial g(u): ui = and < u% < 
1, and u performs periodic oscillatory motion between 
them. When a ^ 0, g(u) generally has three roots {u\ < 
U2 < U3), one of them is necessarily zero. For weak 
fields, the middle root (u 2 ) stays at zero and the molecule 
performs distorted full rotations. When increasing the 
field, a bifurcation happens with the roots of g(u): the 
smallest root u\ becomes stuck at u — 0, and the system 
oscillates in the region < u-i < M3 < 1 where g(u) 
is positive. This corresponds to the so-called pendular 
motion [5[ , when the molecular angular motion is trapped 
by the external field. 

Since molecules experience a time-varying amplitude 
of the optical field while propagating through the deflect- 
ing beam, the total rotational energy of the system and 
the position of the roots ^1,2,3 are changing with time. 
However, these changes are adiabatic with respect to the 
rotational motion, and therefore we can use adiabatic in- 
variants to determine the energy of the system [2l| - |23| . 
The adiabatic invariant related with the coordinate 8 is: 



h = f Pede. (47) 
It is easy to derive from Eas.([2"0]). |03]) and 07]) that : 

(48) 



" 3 VffM 



u(l — u) 



du. 



The energy H of the molecule inside the deflecting field 
as a function of the initial energy H (before entering the 
field) is obtained numerically by solving the Eq.: 



r 



(49) 



where Ig is calculated for a = 0, i.e. in the absence of 
the external field. 

Once the energy of the system H and the polynomial 
g{u) have been found, the average alignment factor is 
simply given by: 



(cos 



J^udu/y/gju) 
Ju'du/^/gJu) 



(50) 



To illustrate the performance of the procedure at real 
experimental conditions [3, 0| , we consider the deflection 
of CS2 molecules at T = 5K (see Fig. [1]), and plot 
the distribution of A at the peak of the deflecting field. 
The results are given in Figs. [TOk and b, for weak (3 • 
10 9 W/cm 2 ) and strong (9- 10 11 W/cm 2 ) deflecting fields, 
respectively. 




FIG. 10: Distribution of A at the peak of the deflecting 
field. Two cases are considered: (a) weak deflecting field 
3-W 9 W/cm 2 and (b) strong deflecting field 9- 10 11 W/cm 2 . In 
(a) we observe the unimodal-rainbow distribution discussed in 
the previous sections. In (b), the effect of the alignment by 
the deflecting field is evident. 



In the case of weak field, we get a unimodular 
rainbow distribution similar to that derived by the 
various methods in the previous sections. In the 
case of strong field we obtain a rotationally-trapped 
distribution, corresponding to the pendular-like mo- 
tion of the molecules at the top of the deflecting pulse [4( . 

To study the deflection of CS2 molecules by a focused 
laser beam, we integrate numerically Eq.© to find the 
deflection velocity. In the integrand of Eq. ([3j) , we substi- 
tute the value of (cos 2 #) calculated by Eq.(f50"l) in every 
point of time. As in the previous sections, we assume that 
x sa v x t (v x = 500to/s) and consider z as a fixed impact 



parameter (z = — 4/im). These assumptions are valid 
even for strong deflecting fields (that align the molecules) 
since the deflection angle is still small. We consider both 
weak and strong deflecting fields, as in Fig. [TQ1 and use 
the values of ojq = 7/xm and r = 14ns (Eq. []} in the 
calculation of the trajectories. 

The distribution of deflected velocities for a thermal 
molecular ensemble (without prealignment) is shown in 
Fig. HU 



40 
I 20 

.Q 



(a) 



0.04 0.08 0.12 0.16 

(b) 
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FIG. 11: Distribution of velocities (or deflecting angles) cal- 
culated from the trajectories of the molecules subject to weak 
(a) and strong (b) deflecting fields. The fields characteristics 
are given in the text. 

In Fig. ITTa (weak field) we essentially verify our as- 
sumption from the previous sections, that the deflection 
in weak fields is linear with A (Eq. QJ. This is seen by 
observing that Fig. [TTb may be indeed obtained by a 
linear transformation of the distribution f(A) from Fig. 
[TOk In Fig. [TTb (strong field) the distribution of deflec- 
tion angles (or of deflection velocities) is still quite broad. 
To our opinion, this results from two different regimes of 
scattering that the molecules experience while traversing 
the deflecting beam: weak deflection at the periphery of 
the beam, and deflection under partial alignment of the 
molecular ensemble in the center of the beam. 

Finally we consider scattering of molecules prealigned 
in the z direction with the pulses having kick strength 
of P = 25. The results are given in Figs. H"2"b and b for 
weak and strong deflecting fields, respectively. 

In the case of weak deflection (Fig. IT2"k). the nar- 
row peak is observed whose nature was already explained 
above. More interestingly, the distribution of the deflec- 
tion angles regained the narrow shape even in the case 
of strong deflection field (Fig. IT2b ) as a result of pre- 
alignment! In our example, the prealignment pulse was 
strong enough to overcome the rotational trapping by the 
deflecting field. As a result, all the molecules performed 
full rotations (but not a pendular motion) despite the 
presence of the strong deflecting field, and we obtain a 
narrow distribution as well. 



o 



V) 

b 

03 
.Q 
O 



300 

200 

100 


0.00 

0.4- 
0.2- 
0.0 




(a) 



0.04 0.08 



0.12 








(b) 



20 
V [m/s] 



40 



FIG. 12: Distribution of deflection velocities for prealigned 
molecules (P — 25). Case (a) is for weak deflecting field, 
and case (b) is for strong deflecting field. In both cases the 
molecules were prealigned parallel to the deflecting field, and 
this prealignment was strong enough to ignore any aligning 
effect of the deflecting field (thus reducing the broadening of 
the deflection, as was explained in the previous sections). 



VI. DISCUSSION AND CONCLUSIONS 

Our results indicate that prealignment provides an 
effective tool for controlling the deflection of rotating 
molecules, and it may be used for increasing the bright- 
ness of the scattered molecular beam. This increase was 
shown both for weak and strong deflecting fields. This 
might be important for nano-fabrication schemes based 
on the molecular optics approach |6j. Moreover, molecu- 
lar deflection by non-resonant optical dipole force is con- 
sidered as a promising route to separation of molecular 
mixtures (for a recent review, see [13] )• Narrowing the 
distribution of the scattering angles may substantially 
increase the efficiency of separation of multi-component 
beams, especially when the prealignment is applied selec- 
tively to certain molecular species, such as isotopes [25| . 
or nuclear spin isomers [26l |27| . More complicated tech- 
niques for pre-shaping the molecular angular distribution 
may be considered, such as confining molecular rotation 
to a certain plane by using the " optical molecular cen- 
trifuge" approach 28j, double-pulse ignited "molecular 
propeller" [29l433| . or planar alignment by perpendicu- 
larly polarized laser pulses [34] - In this case, a narrow 
angular peak is expected in molecular scattering, whose 
position is controllable by inclination of the plane of ro- 
tation with respect to the deflecting field [35| . Laser pre- 
alignment may be used to manipulate molecular deflec- 
tion by inhomogeneous static fields as well (36| (for re- 
cent exciting experiments on post- alignment of molecules 
scattered by static electric fields see [13|)- In particu- 
lar, one may affect molecular motion in relatively weak 
fields that are insufficient to modify rotational states by 
themselves. Moreover, the same mechanisms may prove 
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efficient for controlling inelastic molecular scattering off 
metalic/dielectric surfaces. These and other aspects of 
the present problem are subjects of an ongoing investi- 
gation. 

This research is made possible in part by the historic 



generosity of the Harold Perlman Family. IA is an incum- 
bent of the Patricia Elman Bildner Professorial Chair. 
We thank Yuri Khodorkovsky for providing us the results 
of his related time-dependent quantum calculations. 



[1] H. Stapelfeldt, H. Sakai, E. Constant and P. B. Corkum, 
Phys. Rev. Lett. 79, 2787 (1997); H. Sakai, A. Tarase- 
vitch, J. Danilov, H. Stapelfeldt, R. W. Yip, C. Ellert, 
E. Constant and P. B. Corkum, Phys. Rev. A, 57, 2794 

(1998) . 

[2] B. S. Zhao, H. S. Chung, K. Cho, S. H. Lee, S. Hwang, 

J. Yu, Y. H. Ahn, J. Y. Sohn, D. S. Kim, W. K. Kang, 

and D. S. Chung, Phys. Rev. Lett. 85, 2705 (2000); H. 

S. Chung, B. S. Zhao, S. H. Lee, S. Hwang, K. Cho, S. 

H. Shim, S. M. Lim, W. K. Kang and D. S. Chung, J. 

Chem. Phys. 114, 8293 (2001). 
[3] B. S. Zhao, S. H. Lee, H. S. Chung, S. Hwang, W. K. 

Kang, B. Friedrich and D. S. Chung, J. Chem. Phys. 

119, 8905 (2003). 
[4] S. M. Purcell and P.F. Barker, Phys. Rev. Lett. 103, 

153001 (2009). 

[5] B. Friedrich and D. Herschbach, Phys. Rev. Lett. 74, 
4623 (1995); J. Chem. Phys. Ill, 6157 (1999). 

[6] T. Seideman, J. Chem. Phys. 106, 2881 (1997); J. Chem. 
Phys. 107, 10420 (1997); J. Chem. Phys. Ill, 4397 

(1999) . 

[7] B. Friedrich, Phys. Rev. A 61, 025403 (2000). 

[8] R. J. Gordon, L. Zhu, W. A. Schroeder, and T. Seideman, 

J. Appl. Phys. 94, 669 (2003) 
[9] R. Fulton, A. I. Bishop and P. F. Barker, Phys. Rev. 
Lett. 93, 243004 (2004). 
[10] R. Fulton, A. I. Bishop, M. N. Schneider, P. F. Barker, 

Nature Phys. 2, 465 (2006). 
[11] B. A. Zon and B. G. Katsnelson, Zh. Eksp. Teor. Fiz. 69, 

1166 (1975) [Sov. Phys. JETP 42, 595 (1975)]. 
[12] H. Stapelfeldt and T. Seideman, Rev. Mod. Phys. 75, 
543 (2003). 

[13] V. Kumarappan, S. S. Viftrup, L. Holmegaard, C. Z. 

Bisgaard and H. Stapelfeldt, Phys. Scr. 76 C63 (2007). 
[14] J. J. Larsen, K. Hald, N. Bjerre and H. Stapelfeldt, Phys. 

Rev. Lett. 85, 2470 (2000). 
[15] J.G. Underwood, B. J. Sussman and A. Stolow, Phys. 

Rev. Lett. 94, 143002 (2005); K. F. Lee, D. M. Vil- 

leneuve, P. B. Corkum, A. Stolow and J. G. Underwood, 

Phys. Rev. Lett. 97, 173001 (2006). 
[16] S. S. Viftrup, V. Kumarappan, S. Trippe and H. 

Stapelfeldt, Phys. Rev. Lett., 99, 143602 (2007). 
[17] E. Gershnabel and I. Sh. Averbukh, Phys. Rev. Lett., 

104, 153001 (2010). 
[18] E. Gershnabel, I. Sh. Averbukh and R. J. Gordon, Phys. 

Rev. A 74, 053414 (2006). 
[19] G. B. Arfken, H. J. Weber, Mathematical Methods 

for Physicists, 6th ed. (Elsevier Academic Press, USA, 

2005). 



[20] J. P. Cryan, P. H. Bucksbaum, and R. N. Coffee, Phys. 

Rev. A 80, 063412 (2009). 
[21] H. Goldstein, C. Poole and J. Safko, Classical 

Mechanics, 3rd ed. (Addison Wesley, USA, 2001). 
[22] L. D. Landau and E. M. Lifshitz, Mechanics, 3rd ed. 

(Butterworth-Heinemann, UK, 1976). 
[23] P. Dugourd, I. Compagnon, F. Lepine, R. Antoine, D. 

Rayane and M. Broyer, Chem. Phys. Lett. 336, 511 

(2001). 

[24] B. S. Zhao, Y. M. Koo , D. S. Chung, Analytica Chimica 
Acta 556, 97, (2006). 

[25] S. Fleischer, I. Sh. Averbukh and Y. Prior, Phys. Rev. A 
74, 041403(R) (2006). 

[26] M. Renard, E. Hertz, B. Lavorel, and O. Faucher, Phys. 
Rev. A 69, 043401 (2004). 

[27] S. Fleischer, I. Sh. Averbukh, and Y. Prior, Phys. Rev. 
Lett., 99, 093002 (2007); E. Gershnabel and I. Sh. Aver- 
bukh, Phys. Rev. A, 78, 063416 (2008). 

[28] J. Karczmarek, J. Wright, P. Corkum and M. Ivanov, 
Phys. Rev. Lett. 82, 3420 (1999); D. M. Villeneuve, S. 
A. Aseyev, P. Dietrich, M. Spanner, M. Yu. Ivanov and 
P. B. Corkum, Phys. Rev. Lett. 85, 542 (2000). 

[29] S. Fleischer, I. Sh. Averbukh, and Y. Prior, 
in Ultrafast Phenomena XVI, Proceedings of 
the 16th international conference, p. 75-77, June 
2008, Stressa, Italy, (Springer, Berlin, 2008) ; 
http://www.ultraphenomena.org/programme/tue3.pdf 

[30] S. Fleischer, Y. Khodorkovsky, I. Sh. Averbukh, 
and Y. Prior, in Conference on Lasers and Electro- 
Optics/International Quantum Electronics Conference, 
(Baltimore, Maryland, May 2009) OSA Technical Digest 
(CD) (Optical Society of America, 2009), paper ITuC5. 

[31] S. Fleischer, Y. Khodorkovsky, Y. Prior, and I. Sh. Aver- 
bukh, New J. Phys. 11, 105039 (2009). 

[32] A.G. York, Opt. Express, 17, 13671 (2009). 

[33] Kenta Kitano, Hirokazu Hasegawa, and Yasuhiro 
Ohshima, Phys. Rev. Lett., 103, 223002 (2009). 

[34] M. Lapert, E. Hertz, S. Guerin and D. Sugny, Phys. Rev. 
A 80, 051403(R) (2009). 

[35] J. Floss, E. Gershnabel and I. Sh. Averbukh, to be pub- 
lished 

[36] E. Gershnabel and I.Sh. Averbukh, to be published 

[37] L. Holmegaard, J. H. Nielsen, I. Nevo and H. Stapelfeldt, 

Phys. Rev. Lett. 102, 023001 (2009); F. Filsinger, J. 

Kiipper, G. Meijer, L. Holmegaard, J. H. Nielsen, I. Nevo, 

J. L. Hansen and H. Stapelfeldt, J. Chem. Phys. 131, 

064309 (2009). 



